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Abstract. Empirical fitting functions, describing the be- 
haviour of the A4000 A break, D4000, m terms of effec- 
tive temperature, metallicity and surface gravity, are pre- 
sented. For this purpose, the break has been measured in 
392 stars from the Lick/IDS Library. We have followed a 
very detailed error treatment in the reduction and fitting 
procedures, allowing for a reliable estimation of the break 
uncertainties. This calibration can be easily incorporated 
into stellar population models to provide accurate predic- 
tions of the break amplitude for, relatively old, composite 
systems. 
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1. Introduction 

The spectroscopic study of the blue and near-UV region 
around A4000 A has proven to be a useful tool to investi- 
gate the stellar populations of composite stellar systems. 
Obviously, this spectral region is specially suited to detect 
the presence of young stars and therefore, to study star 
formation histories. Furthermore, in a pioneering work, 
Morgan (1959) already showed that the difference in in- 
tensity of the continuum level on the two sides of H£ 
(A3889A) in galactic globular clusters correlated well with 
metal abundance. 

Although different absorption line-strength indices 
have been defined in this spectral range to understand the 
stellar composition of early-type objects (e.g., Faber 1973; 
Burstein et al. 1984; Rose 1984, 1994; Tripicco 1989; Wor- 
they et al. 1994, hereafter W94; Jones & Worthey 1995; 
Worthey & Ottaviani 1997; Vazdekis & Arimoto 1999), 
some of them are quite dependent on spectral resolution 
(and therefore velocity dispersion) and, in many cases, 
their use requires relatively high signal-to-noise ratios. In 
this sense, an interesting spectral index that avoids these 
problems is the A4000 A break. We have already demon- 
strated (Cardiel et al. 1998) that this discontinuity can be 
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measured with a relative error of ~ 10% with a signal-to- 
noise ratio per A ~ 1. Thus, the break is well suited to be 
measured in faint objects or at low surface brightnesses. 
However, this advantage, due to the large wavelength in- 
terval employed in its definition, also translates into an 
important drawback: many absorption lines are included 
in the break bandpasses. Therefore, the behaviour of the 
break is expected to be complex. 

In this work we are using the definition adopted by 
Bruzual (1983), who defined this spectral index as the 
ratio of the average flux density, F v (erg s -1 cm~ 2 Hz -1 ), 
in two bands at the long- and short-wavelength side of 
the A4000 A discontinuity, in particular 
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where (Af, A J, Xf, A J) = (3750, 3950, 4050, 4250) A. Ex- 
cept for the combination of v and A, (due to the measure- 
ment method employed by Bruzual, in which the break 
was obtained from galaxy spectra plotted as F v versus A), 
and for not being in logarithmic units, this definition re- 
sembles that of a color. Because of this, the D4000 can be 
considered as a pseudo-color. However, it must be pointed 
out that, compared to the classical color indices in this 
spectral range (like U — B), the 4000 A break is much 
less sensitive to extinction by dust, and hence, it is more 
valuable to investigate the stellar content of galaxies. In 
particular, using the average interstellar extinction curve 
from Savage & Mathis (1979), the D4000 can be corrected 
from internal reddening (and galactic extinction for ob- 
jects at zero redshift) using the following expression 
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It is interesting to note that the break definition is still 
the same for redshifted objects, where the integral limits 
must be properly modified, and that absolute fluxes are 
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not required. As a working expression, the 4000 A break 
can be rewritten as 
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where fx is the flux density per unit wavelength (in ar- 
bitrary units), and z is the redshift of the object being 
measured. 

Up to date, many authors have employed the D4000 to 
study the stellar composition and star formation history of 
early-type galaxies (e.g. McClure & van den Bergh 1968; 
Spinrad 1980, 1986; Bruzual 1983; Laurikainen & Jaakkola 
1985; Hamilton 1985; Dressier 1987; Dressier & Shectman 
1987; Johnstone et al. 1987; Kimble et al. 1989; Dressier & 
Gunn 1990; Rakos et al. 1991; Chariot et al. 1993; Chariot 
& Silk 1994; Davidge & Clark 1994; Songaila et al. 1994; 
Bclloni et al. 1995; Davidge & Grinder 1995; Cardiel et al. 
1995, 1998; Abraham et al. 1996; Hammer et al. 1997; Bar- 
bara & Poggianti 1997; Longhetti et al. 1998; Ponder et 
al. 1998). The reliable analysis of the break measurements 
rests on the comparison of the data with the predictions 
of stellar population models (e.g. Worthey 1994; Bruzual 
& Chariot 1996). So far, such predictions are computed by 
using either model atmospheres, or stellar libraries with a 
poor coverage of the atmospheric parameter space, espe- 
cially in metallicity. 

In this paper we present an empirical calibration of 
the A4000 A break as a function of the main atmospheric 
stellar parameters (namely effective temperature, surface 
gravity and metallicity) in an ample stellar library which 
covers an appropriate range of parameters to study rela- 
tively old stellar populations. One of the main advantages 
of using fitting functions to describe the behaviour of spec- 
tral indices is that they allow stellar population models to 
include the contribution of all the required stars, through 
a smooth interpolation in the space defined by the fitted 
stellar parameters. The usefulness of this approach has 
been demonstrated by the successful inclusion of similar 
fitting functions in recent evolutionary synthesis models 
(eg. Worthey 1994; Vazdekis et al. 1996; Bressan et al. 
1996; Bruzual & Chariot 1996). 

It is important to keep in mind that the empirical cali- 
bration is only a mathematical representation of the break 
behaviour as a function of atmospheric stellar parameters, 
and that we do not attempt to obtain any physical justi- 
fication of the derived coefficients. 

We briefly review the previous works devoted to un- 
derstand the D4000 in section 2. The star sample is given 
in section 3. The observations and data reduction are de- 
scribed in section 4, whereas section 5 contains a descrip- 
tion of the error analysis. In section 6 we show the be- 
haviour of the measured D4000 values as a function of the 
stellar atmospheric parameters. The fitting strategy and 



the resulting empirical function are presented in sections 7. 
Finally, in section 8 we give a summary, providing a public 
FORTRAN subroutine written by the authors to facilitate 
the computation of the D4000 using the fitting function 
presented in this paper. Sections 4 and 5 are rather tech- 
nical, due to the inclusion of a lengthy explanation of the 
data and error handling. We suggest the reader not in- 
terested in such details to scan Tables |l| and ^, and skip 
those sections. 



2. Previous works: understanding the D4000 

The A4000 A break is a sudden onset of absorption fea- 
tures bluewards 4000 A which is clearly noticeable for 
stellar types cooler than GO (see Fig. ||). In Fig. [l] we 
show a typical spectrum of a cool star in this spectral 
region, together with the identification of the most promi- 
nent spectral features. Considering the large wavelength 
range employed in the measurement of the D4000 , it is ex- 
pected the strength of this discontinuity to be a function 
of the distribution of the continuum light in this region 
(governed by the effective temperature) modulated by the 
absorption line strengths (which must depend primarily on 
both temperature and metallicity, and secondly on grav- 
ity). This behaviour converts the break in a potential tool 
to investigate composite stellar populations in early-type 
systems. 

The relevance of the line-blanketing discontinuity near 
A4000 A was the object of a systematic study by Wildey et 
al. (1962). These authors measured the energy subtracted 
in the spectra of some stars due to Fraunhofer lines, show- 
ing that the effect was important below A4000 A. Van den 
Bergh (1963), and van den Bergh & Sackmann (1965) de- 
fined a break, A, as the ratio of the smoothed observed 
continuum at both sides of A4000 A. These authors mea- 
sured this break in a sample of 200 stars concluding that 
A depended both on stellar metallicity and B— V color. 
Analogous discontinuity definitions, like C(38 — 41) (Mc- 
Clure & van den Bergh 1968), and T(38 - 41) (Carbon et 
al. 1982), have been also employed in the spectroscopic 
analysis of stars, star clusters and galaxies. 

Using spectrophotometric stellar libraries, Bruzual 
(1983) and Hamilton (1985) studied the variation of the 
A4000 A break with spectral types and luminosity classes 
(compare Fig. 3 in Bruzual with Fig. 6 in Hamilton). 
Both authors concluded that, as a function of tempera- 
ture, the D4000 increases slowly for spectral types in the 
range from 05 to GO, and faster from GO to M0, whereas 
the break decreases for the later types, M0 to M5. In ad- 
dition, whilst for spectral types hotter than GO the break 
does not depend on gravity, a clear dichotomy between 
main sequence stars on one hand, and giant and super- 
giant stars on the other, is apparent for lower tempera- 
tures. Given the scarcity of the employed stellar libraries, 
no dependence on metallicity could be obtained in these 
works. 
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Fig. 1. Spectrum of the star HD72324 (G9 III) in the re- 
gion around the A4000 A break. As a reference, we have 
also plotted the most intense (EW > 200 mA) Fraunhofer 
lines from the sun (using the tabulated data from Lang 
1974 — original source Moore et al. 1966 — ), together with 
the Balmer lines and two CN molecular bands (the central 
bandpass of the CN3883 index defined by Pickles 1985, 
and the absorption band of the S(4142) index employed 
by Smith et al. 1997) which can be found in this spec- 
tral range. The contribution of the atomic metallic lines, 
especially from Fe I and Mg I, becomes very important 
bluewards A4000 A. 



From the analysis of moderate-resolution spectra of 
950 galaxies in 12 rich clusters, Dressier & Shectman 
(1987) argued that, in composite stellar populations, the 
break is insensitive to changes in metal abundance, at least 
in the metallicity range spanned by their galaxy sample. 
This result was employed by Munn (1992) to conclude that 
the CN-D4000 diagram is effective at separating metallic- 
ity and age effects on the integrated spectra of early-type 
galaxies. However, Kimble et al. (1989) obtained that the 
break correlated strongly with metallicity indicators, such 
as the Mgb index. 

More recently, Poggianti & Barbaro (1997), working 
with Kurucz's models, have obtained a theoretical cali- 
bration of the break as a function of stellar parameters. 
They present (Fig. 1 in their paper) the behaviour of the 
D4000 in the ranges 5500 < T eS < 35000 K, < \ogg < 5, 
and — 2 < \ogZ < 0. This work clearly shows the strong- 
dependence of the break on effective temperature, as pre- 
viously reported from the studies based on stellar libraries, 
and quantifies, for a small sample of temperatures, the 
variation of the break as a function of metallicity and grav- 
ity. The D4000 is shown to be insensitive to metallicity for 
hot stars (T c g — 9000 K), whereas the contrary is true 
for T ff = 5500 K. In addition, using the stellar spectra of 
Straizys & Sviderskiene (1972; note that these spectra are 
those also employed by Bruzual 1983), these authors ob- 



tain that, for stars with 3500 < T c g < 5500 K, the D 40 oo 
always exhibits values above 2, with a maximum of 3 at 
T e ff = 4000 K. Using this theoretical calibration, Barbaro 
& Poggianti (1997) have also elaborated an evolutionary 
synthesis model which predicts, in the integrated spec- 
trum of a galaxy, the variation of the D4000 as a function 
of the star formation rate (SFR). More interestingly, they 
conclude that the break can be employed to yield the ra- 
tio of the SFR averaged over the last 5 billion years to the 
present SFR. 

From all these previous works, it is quite clear that the 
D4000 is a suitable tool to face the study of stellar systems, 
in particular to reveal their stellar composition. However, 
a detailed empirical calibration, such as that presented 
in this paper, is needed to i) overcome the unavoidable 
uncertainties associated to the theoretical calibrations, ii) 
extend our understanding of the break behaviour for stars 
with T e ff < 5500 K (note that these late-type stars con- 
stitute a fundamental ingredient in the modeling of old 
stellar populations), and iii) use in conjunction with other 
indices previously calibrated with the same stellar library. 

3. Star sample 

In order to derive a confident empirical calibration of the 
D4000, we decided to measure this spectral feature in all 
the stars belonging to the Lick/IDS Library (Burstein et 
al. 1984; Faber et al. 1985; Burstein et al. 1986; Gorgas 
et al. 1993, hereafter G93; and Worthey et al. 1994). The 
suitability of the Lick/IDS library to obtain empirical fit- 
ting functions has been widely demonstrated by the works 
of G93, W94, and Worthey & Ottaviani (1997), who, in 
overall, have derived analytical expressions for 25 spectral 
indices in the 4000-6000 A region. Trager et al. (1998), 
and references therein, have extensively shown the useful- 
ness of the Lick/IDS absorption- line index system in the 
study of old stellar populations. In brief, the Lick/IDS li- 
brary contains 460 stars of all spectral types and luminos- 
ity classes. Although a large fraction are field stars from 
the solar neighbourhood, members of open clusters (cov- 
ering a wide range of ages) and galactic globular clusters 
(with different metallicities) are also included. 

Table [l] lists the final data sample, which comprises 392 
stars out of the original set of 424 objects with published 
atmospheric parameters in W94 (the remaining 32 stars 
presented high uncertainties in the derived measurements 
and were not included in order to guarantee the quality 
of the final dataset) . The first two columns list the Henry 
Draper Catalog number, if any, and other common desig- 
nation. The spectral type and luminosity class are given in 
the third column (see references at the end of the table). 
The following three columns list the atmospheric parame- 
ters employed in this work. Most of them are those listed 
in the electronic table of W94. Since we wish to keep the 
consistency with the previous fitting functions for the Lick 
indices, we have been reluctant to change any atmospheric 
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Fig. 2. Gravity-temperature diagram for the sample of 
stars used to derive the empirical fitting function. Differ- 
ent symbols are used to indicate stars of different metal- 
licities, as shown in the key. 



parameter. However, we have introduced new parameters 
for some of the stars which lack them in W94, and, in some 
cases, the original parameters have been improved using 
more recent or reliable data from the literature (see the 
description in the table notes). As a whole, effective tem- 
perature, gravity and metallicity data are finally available 
for 383, 382 and 354 stars, respectively. 

To summarize, the Lick/IDS library spans the follow- 
ing ranges in atmospheric parameters: T e g from 2747 to 
21860 K, log. 9 from to 5.12, and [Fe/H] from -2.70 to 
0.46. It is important to note that the star sample does 
not homogeneously cover the parameter space, being more 
densely populated in the metallicity interval —1.0 dex to 
+0.5 dex, whereas the gravity range is wider close to solar 
metallicity. A more important deficiency of the library is 
the paucity of hot stars. Since we are also interested in 
applying the derived fitting functions to the modeling of 
younger stellar populations, we have enlarged the sample 
with the inclusion of additional hot stars. For this purpose, 
we have employed flux-calibrated spectra corresponding to 
representative spectral types and luminosity classes from 
the compilation of Pickles (1998). In particular, 43 spec- 
tra of 05-FO stars with 6720 < T cS < 39800 K have 
been incorporated. It is worth noting that we have checked 
that no systematic offset is apparent between our D4000 
measurements and those in Pickles' spectra. To display 
the range of stellar parameters covered by the resulting 
library, in Fig || we present a log g-T e g diagram of the 
calibrating stars. 



4. Observations and data reduction 

The stellar library was observed in a total of 13 observing 
runs from 1991 to 1996 using the JKT, INT and WHT 
at the Roque de los Muchachos Observatory (La Palma, 
Spain), the 3.5 m telescope at Calar Alto Observatory 
(Almcria, Spain), and the 2.12 m telescope at San Pe- 
dro Martir Observatory (Mexico). The bulk of the spectra 
(320 out of a total of 650 spectra) were taken in run 6, 
where the Richardson-Brealey Spectrograph was used to 
obtain 2.8 A resolution spectra in the blue spectral range. 
A description of each observing run, including relevant ob- 
servational configuration parameters, is given in Table § 
The number of stars observed in each run is quite variable 
since some of the runs were not devoted to the calibra- 
tion of the break and, in these cases, only bright stars of 
the library were observed during twilight periods. The last 
column in Table [l] list the run numbers in which each star 
was observed. 

The reduction of the data was performed with our 
own reduction package REPmEQ (Cardiel & Gorgas 1999), 
which allows a parallel treatment of data and error spectra 
(see below). We followed a standard reduction procedure 
for spectroscopic data: bias and dark subtraction, cos- 
mic ray cleaning, flat-fielding, wavelength calibration, C- 
distortion correction, sky subtraction, atmospheric extinc- 
tion correction and flux calibration (we did not attempt to 
obtain absolute fluxes since, as most line-strength indices, 
the break only requires relative fluxes) . Cluster stars spec- 
tra were also corrected from interstellar reddening, using 
the color excesses quoted in Tables 4 and A3 from G93 
and W94, respectively, and the averaged extinction curve 
of Savage & Mathis (1979). 

In order to optimize the observing time during most 
of the runs, comparison arc frames were not taken next to 
each star observation. Instead of this, we only acquired arc 
exposures for a selected subsample of stars, which com- 
prised objects with a complete coverage of all the spec- 
tral types and luminosity classes observed in each run. 
The wavelength calibration of the rest of the stars was 
performed by comparing them with the reference spectra. 
The repetition of this procedure with different reference 
spectra allow us to guarantee that wavelength calibration 
errors were always ^ 0.2 pixels. 

5. Random errors and systematic effects 

Since the aim of this paper is to derive an analytical rep- 
resentation of the behaviour of the D4000 as a function 
of effective temperature, metallicity and surface gravity, 
the sources of error are two-fold. In one hand, an impor- 
tant error source are the uncertainties in the adopted at- 
mospheric stellar parameters. Detailed discussions of the 
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Table 1. Data sample, stellar atmospheric parameters employed to derive the fitting function, and break measure- 
ments. 
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4368 


1.200 


-1 


.50 


2 


.651 


0.084 


0.010 


6(^) 




Ml 3 A 1 71 

1V1 1 A 111 


AGB 


4401 


1.070 


-1 


.50 


3 


.159 


0.071 


0.577 


6(^) 




Ml 3 R 786 

1V11U JJ I LJU 


GB 


3963 


0.620 


-1 


,50 


3 


.401 


0.101 


—0.377 


6(1) 




Ml 3 R 81 8 

1V1 1 <J JJ Ulu 


AGB 


5089 


1.890 


-1 


.50 


I 


.691 


0.089 


—0.180 


7 




TV/TO QQO 


Lz-D 


A Pil 9 




-1 


.70 


o 

z 


3S9 


n 1 1 i 
U.lll 


U.UDU 


7 




M3 111-28 


GB 


4122 


0.730 


-1 


.70 


2 


.708 


0.072 


-0.436 


6(4) 




M3 IV-25 


GB 


4355 


1.210 


-1 


.70 


2 


.516 


0.122 


-0.062 


7 




M5 1-45 


HB 


5850 


2.610 


-1 


.30 


1 


.518 


0.087 


-0.095 


7 




M5 11-51 


GB 


4569 


1.690 


-1 


.30 


2 


.444 


0.117 


0.021 


7 




M5 11-53 


HB 


10460 


3.660 


-1. 


.30 


1 


.347 


0.077 


0.046 


7 




M5 11-76 


HB 


6060 


2.700 


-1 


.30 


1 


.422 


0.060 


-0.136 


7(2) 




M5 111-03 


GB 


4095 


0.660 


-1 


.30 


3 


.626 


0.127 


0.258 


6(2) 




M5 IV- 19 


GB 


4139 


0.840 


-1 


.30 


3 


.204 


0.102 


-0.039 


6(2) 




M5 IV-59 


GB 


4300 


0.850 


-1 


.30 


3 


.056 


0.105 


0.193 


6(2) 




M5 IV-86 


HB 


5429 


2.460 


-1 


.30 


1 


.509 


0.060 


-0.253 


7(2) 




M5 IV-87 


HB 


5950 


2.630 


-1 


.30 


1 


.572 


0.087 


-0.014 


7 



6 
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Table 1. Continued. 



Name 


Type" 


rji b 

J-cB 


log/ 


[Fe/H] b 


D4000 


AD4ooo c 


Res. d 


Runs' 




HR 


i± i Do 


9 ^90 


o i n 

— U. 1U 


9 Q^P 
Z.yoD 


n 0Q8 
u.uyo 


081 
U.Uol 


D 


Mf\7 F 0Q4 


O '. ' / ) 


(SI m 

U1UO 


4 070 


—0 i n 

U. 1U 


1 (SQ8 


n 0Q9 


078 
u.u 1 


p 
u 


IVIU I -T 1UO 


CR 




9 9^0 


o in 

U. 1U 


^ 41 1 
0.111 


n 1 09 

U. lUi 


1 ^8 

U. IOO 


fi 
u 


\I\(\7 P 1 08 
IVIO ( r lUo 


CR 

Kill 




l.OOU 


o i n 

— U. 1U 


789 
O. ( OZ 


U. 1UO 


91 R 
U.Z1D 


p 


A/TRY piic; 

IVIU I -T _L _LO 


^CR 




8Q0 


o i n 

U. 1U 


1 7P^ 
1 . 1 uo 


n no*} 

u.uyo 


nfifi 

u.uuu 


u 


A/TRT 1? 117 

IVIU l -T ± ± i 


ICR 




7Q0 
o. / yu 


o in 

U. 1U 


1 80S 


n Ofifi 
u.uuu 


491 
— u.iz 1 




MP\7 F 1 1 Q 
iviu f -T _l _i_y 


ij Ur-D 


UUOl 


Q40 


o in 

— U. 1U 


1 fi83 
1 .uoo 


n 0Q9 
u.uyz 


n9fi 

u.uzu 


p 
u 


MR? F 1 9^ 

IVIU f -T J_^<J 


ij Ur-D 


uuuu 


4 "^40 


o i n 

U. 1U 


1 P^KP 
1 .uou 


n 0^7 

u.uo / 


nni 

u.uuo 


0(6) 


IVIU t C 1U4: 


HB 


4707 

i± 1 KJ 1 


2.220 


—0 i n 


^ ni ^ 

O.UIO 


n 0Q8 

u.uyo 


nss 

u.uoo 


p 
u 


IVIU I -T _L 1 U 


CR 

\ 1 1 1 


A9QA 


1 S^O 

l.OOU 


o in 

U. 1U 


^ 09^ 
o.yzo 


n 1 0^ 

U. 1UO 


41 4 


« 
u 


A/TRy P 1 7^ 
iviu / r ± / o 


ij KjJD 


uuuu 


4 "340 


o i n 

U. 1U 


1 RRI 


n 0^7 

u.uo / 


01 ^ 

U.UIO 


0(6) 


iviu i r i yo 


(1R 
\iii 


i±iyOO 


^00 


—0 i n 

U. 1U 


9 7^iR 
z. 1 ou 


n oro 

u.uuu 


1 9^ 

U. liO 


Pd) 

0(6) 


■N/Tfi7 P 99/1 


CRf 


/171 7 


9 Pi^O 


o i n 

— U. 1U 


0. 1U1 


n OQQ 
u.uyy 


1 8Q 

u. ioy 


P 
D 


Mfi7 p OQ1 
IVIU I -T ZrO-L 


CR 
iii) 




9 Q^O 
z. you 


o i n 

— U. 1U 


9 88^ 

i.OOO 


n 071 

u.u t 1 


1 Qfi 
u. iyu 


0(A) 


A/TR7 T 1 7 

IVIU ( 1 J. ( 


tlR 
iii> 


4Q77 


^ ^70 
o.o/u 


—0 in 

U. 1U 


9 7^1 

i. / Ol 


n ofii 

U.UUl 


1 ^l 
u. 101 


0(6) 


A/TRy TT 99 

IVIU 1 11-ZrZr 




^070 

OU 1 u 


o.uou 


n in 

U. 1U 


9 ^41 


n ofifi 

u.uuu 


048 
u.uio 


0(6) 


\/\a*7 TV 90 

IVIU i 1 V - ZrU 


CR 




9 7^0 
Z. ( ou 


o i n 

— U. 1U 


9 Q49 


n 070 

u.u t u 


npz, 

u.uuo 


0(A) 


A/TR7 TV fS8 

IVIU / 1 V uo 


ICR 

1 J I. J i 1 




7^0 
O. 1 ou 


—0 in 

U. 1U 


9 ^nn 


n ofio 

u.uuu 


—0 118 
u. 1 10 


0(6) 


TV/TR7 TV 77 

IVIU I IV-// 


ICR 


OUiO 


^70 
o.o / u 


n in 

U. 1U 


9 ^R7 


n ofifi 
u.uuu 


nno 

u.uuy 


0(6) 


A/TR7 TV 81 

IVIU I 1 V - 0-L 


ICR 


ooou 


1 7R0 
o. / uu 


o i n 

— U. 1U 


9 1 8R 
z. IOU 


n o^fi 

u.uou 


nn^ 

u.uuo 


0(6) 


M71 1 OQ 


AGB 


4(SQ7 


1 (S70 

1.UIU 


—0 z,fi, 


^ ^P 
o.oou 


n 9^1 


R8R 
u.uou 


\\ 


A/T71 1 91 


CR 

iiii 


AA\ A 


1 4R0 


^Pi 
u.ou 


9 874 
z.o / ^1 


n 1 qi 

u. iy 1 


1 R8 

U. 1UO 


1 1 

1 1 


M71 1 *31 
1V1 (1 i-oi 










1 798 
1. / Zo 


n 1 99 

U. 1ZZ 


T 


1 1 

1 1 


A/T71 1 ^4 

IVI f _L _L 04: 


HB 


^0(19 

OKJV)^ 


9 470 


—0 z,P, 
u.ou 


1 Q' I i' : ! 
1 . yoo 


n 1 

u. 100 


— H4fi 


11 


1V1 / 1 i-ou 










9 497 


n 1 ^9 

U. lOi 


T 


1 1 

1 1 


M71 1 "V7 
1V± / ± _l-o 1 


CR 

iii) 




9 1 80 
z. lou 


W 
u.ou 


9 nns 


n 1 c i' : ! 

U. IOO 


7R9 
u. / uz 


1 1 

1 1 


A/T71 1 


HB 




9 4^0 


—0 z,P, 
u.ou 


1 7^!9 
1 . ( 0^ 


0.141 


— ^QS 
u.oyo 


11 


A/T71 1 zLI 

1V1 ( ± _L~4:_L 


HR 

III) 


ouyo 


9 480 


u.ou 


1 84n 


n 1 94 

U. Iii 


408 


1 1 

1 1 


A/T71 1 ^ 


CR 




1 490 


u.ou 


4 9Rn 


n 9Q9 


8^ 
u.000 


1 1 

1 1 


M71 1 


CR 

III) 




2.440 


—0 z,P, 
u.ou 


9 sns 

i.OUO 


n ^^fi 

u.oou 


0.112 


11 


A/T71 1 fS*3 


ACR 
. ) i. i i > 


471 7 


1 890 


u.ou 


9 ^98 


n im 

U. IOU 


noo 

u.uyy 


1 1 

1 1 


A/T71 1 (\A 

IVI i 1 l-O^i 


CR 

Kill 




1 Pil 

I.OIU 


— U.OD 


9 8R^ 
z.oDo 


n 909 




— u.000 


1 1 

1 1 


A/T71 1 fS'S 


CR 


4640 


2.200 


—0 ^ifi 
u.ou 


9 nsn 


n 91 


— fi41 


11 


A/T71 1 fSfS 


ACR 

. ) K 1 1 ) 


^yu 


1 470 

1 .4: / U 


^Fi 
u.ou 


4Q9 
o.iyz 


u.000 


^7^ 
u.o 1 


1 1 

1 1 


A/T71 1 71 
IVI ( ± _l- i ± 


CR 




1 81 

I.OIU 


^Pi 
u.ou 


9 ^4n 


n 1 Q7 

u. iy / 


P^,7 
u.uo I 


1 1 

1 1 


A/T71 1 7^ 


CR 

Kill 


4846 


2.520 


—0 ^P, 
u.ou 


9 1 ^R 

Z. IOU 


n 1 fii 

U. 1U1 


— ^97 

U.O^ ( 


11 


A/T71 1 Qt\ 
ivi / ± i-yo 


ACR 

. 1 Kit) 




1 R70 
1 .u / u 


^P. 
u.ou 


^117 
0. 1 1 1 


n 1 80 

U. IOU 


All. 


1 1 

1 1 


A/T71 1107 
IVI / 1 _lu / 


ACR 


4Q77 


1 Q^O 
i . you 


— u.ou 


9 ^87 


U. 1UO 


u.zoo 


1 1 

1 1 


A/T71 1 1 0Q 
ivi ( i _i_ iuy 


CR 
it I ) 




9 ^70 


—0 ^P, 
u.ou 


9 71 Q 
i. / iy 


n ^00 

u.ouu 


0.142 


11 


A/T71 A 9 

IVI i 1 rtZr 


HR 
ill) 




9 ^70 


^P 
u.ou 


9 Q1 8 


n 41 1 

U.111 


^9Q 
u.ozy 


1 1 

1 1 


A/T71 A A 

IVI / 1 .ty^ 


ACR 
. ) lit) 


4098 


740 

U. 1 4:U 


^P 
u.ou 


^ PV) 
o.uoz 


n 1 84 

u. 101 


1 pz, 

U. 1UO 


1 1 

1 1 


M71 Afi 

IVI ( 1 rtU 










^ R4Pi 


n 801 

U.OU1 


+ 
T 


11 


M71 A 7 

IVI i 1 r\ ( 










Q49 


n ^fip; 

u.ouo 


T 


1 1 

1 1 


M71 A9 


GB 


4156 


1.390 


—0.56 


4.119 


0.430 


0.607 


11 


M71 C 


HB 


4965 


2.390 


-0.56 


2.177 


0.186 


-0.189 


11 


M71 S 


GB 


4252 


1.390 


-0.56 


3.824 


0.352 


0.501 


11 


M71 X 


HB 


5240 


2.490 


-0.56 


2.499 


0.236 


0.362 


11 


M92 I- 10 


HB 


9290 


3.440 


-2.20 


1.327 


0.076 


-0.055 


7 


M92 1-13 


HB 


5680 


2.220 


-2.20 


1.498 


0.083 


0.022 


7 


M92 11-23 


HB 


7510 


3.050 


-2.20 


1.437 


0.086 


-0.046 


7 


M92 111-13 


GB 


4228 


0.580 


-2.20 


2.473 


0.077 


-0.188 


6,7 


M92 IV- 114 


GB 


4728 


1.580 


-2.20 


1.702 


0.090 


-0.188 


7 


M92 VI- 74 


HB 


5950 


2.210 


-2.20 


1.535 


0.085 


0.101 


6(2) 


M92 IX- 12 


AGB 


5477 


1.930 


-2.20 


1.771 


0.057 


0.249 


9,12 
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Table 1. Continued. 



HD 


Name 


Type" 


T c e b 


log/ 


[Fe/H] 6 


D4000 


AD4ooo c 


Res. 


Runs' 




M92 XII-8 


GB 


4488 


1.000 


—2.20 


1.843 


0.069 


—0.316 


6(2) 




M92 XII-24 


HB 


11100 


3.750 


-2.20 


1.498 


0.060 


0.235 


9,12 




NGC188 1-20 










2.605 


0.060 


f 
1 


10 




NGC188 1-55 


SGB 


5440 


3.910 


0.00 


2.076 


0.043 


—0.063 


9,10 




NGC188 1-57 


GB 


4707 


3.070 


0.00 


3.035 


0.047 


0.029 


9,10 




NGC188 1-61 


GB 


4930 


3.350 


0.00 


2.984 


0.061 


0.245 


10 




NGC188 1-69 


GB 


4401 


2.350 


0.00 


3.479 


0.044 


0.066 


9,10 




NGC188 1-75 


GB 


4906 


3.220 


0.00 


2.879 


0.060 


0.115 


10 




NGC188 1-85 


GB 


4949 


3.550 


0.00 


2.820 


0.062 


0.102 


10 




NGC188 1-88 


SGB 


5285 


3.850 


0.00 


2.358 


0.046 


0.059 


9,10 




NGC188 1-97 


SGB 


5213 


3.820 


0.00 


2.577 


0.053 


0.196 


9,10 




NGC188 1-105 


HB 


4585 


2.190 


0.00 


3.221 


0.043 


0.067 


9,10 




NGC188 1-116 


? 


5070 


3.230 


0.00 


2.135 


0.042 


—0.430 


9,10 




NGC188 11-52 


SGB 


5532 


3.940 


0.00 


2.137 


0.061 


0.085 


10 




NGC188 11-64 


SGB 


5837 


4.090 


0.00 


1.696 


0.209 


-0.113 


13 




NGC188 11-67 


SGB 


5935 


4.130 


0.00 


1.726 


0.207 


-0.017 


13 




NGC188 11-69 


SGB 


5972 


4.170 


0.00 


1.712 


0.211 


-0.007 


13 




NGC188 11-72 


GB 


4387 


2.410 


0.00 


3.507 


0.234 


0.074 


13 




NGC188 11-76 


HB 


4553 


2.180 


0.00 


3.276 


0.220 


0.078 


13 




NGC188 11-79 


SGB 


5138 


3.550 


0.00 


2.405 


0.220 


-0.065 


13 




NGC188 11-88 


GB 


4520 


2.710 


0.00 


2.165 


0.219 


-1.078 


13 




NGC188 11-93 


SGB 


5503 


3.930 


0.00 


2.025 


0.219 


-0.054 


13 




NGC188 11-122 


GB 


4949 


3.410 


0.00 


2.788 


0.250 


0.072 


13 




NGC188 11-126 


GB 


4949 


3.450 


0.00 


3.047 


0.269 


0.330 


13 




NGC188 11-181 


GB 


4276 


2.190 


0.00 


3.491 


0.068 


-0.107 


9,13 




NGC188 11-187 


GB 


4949 


3.330 


0.00 


3.007 


0.067 


0.294 


9,13 




NGC6171 04 


HB 


6100 


2.750 


-0.99 


2.915 


0.932 


1.313 


12 




NGC6171 45 


HB 


5920 


2.840 


-0.99 


1.627 


0.649 


-0.030 


12 




NGC7789 338 










2.635 


0.079 


t 


6(2) 




NGC7789 415 


GB 


3859 


1.060 


-0.12 


4.421 


0.126 


0.243 


6(2) 




NGC7789 468 


GB 


4228 


1.600 


-0.12 


3.941 


0.142 


0.339 


6(2) 




NGC7789 489 










3.157 


0.086 


t 
1 


6(2) 




NGC7789 501 


GB 


4085 


1.380 


-0.12 


3.961 


0.177 


0.123 


6(2) 




NGC7789 575 


GB 


4512 


1.980 


-0.12 


2.907 


0.219 


-0.265 


13 




NGC7789 669 


GB 


4216 


1.570 


-0.12 


4.053 


0.135 


0.431 


6(2) 




NGC7789 676 


HB 


4988 


2.320 


-0.12 


2.696 


0.216 


0.108 


13 




NGC7789 859 


GB 


4626 


2.270 


-0.12 


2.517 


0.212 


-0.499 


13 




NGC7789 875 


HB 


4965 


2.360 


-0.12 


2.517 


0.212 


-0.095 


13 




NGC7789 897 


HB 


4965 


2.350 


-0.12 


2.064 


0.208 


-0.548 


13 




NGC7789 971 


GB 


3831 


1.030 


-0.12 


4.403 


0.176 


0.264 


6(2) 


1461 


HR 0072 


GO V 


5932 


4.300 


0.29 5 


2.364 


0.099 


0.519 


5 


2665 




G5 III 


5000 


2.500 


-2.00 


1.518 


0.090 


-0.243 


6 


2857 


BD-06 0086 


A2 


7563 1 


2.670 1 


-1.58 1 


1.629 


0.091 


0.144 


6 


3546 


HR 0163 


G8 III 


4957 


2.600 


-0.66 


2.247 


0.055 


-0.075 


3 


3567 




F5 VI 


6001 


4.160 


-1.44 


1.499 


0.084 


0.120 


5 


3651 


HR 0166 


K0 V 


5487 


4.520 1 


-0.18 1 


2.741 


0.093 


0.749 


5 


4307 


HR 0203 


GO V 


5650 


3.930 


-0.52 


1.864 


0.090 


0.144 


5 


4614 


HR 0219 


GO V 


5749 


4.420 4 


-0.20 4 


2.000 


0.134 


0.211 


4 


4628 


HR 0222 


K2 V 


4960 


4.600 1 


-0.29 1 


2.758 


0.091 


0.209 


5 


4656 


HR 0224 


K5 III 


3953 


1.100 


-0.07 


4.784 


0.113 


0.713 


3 


5395 


HR 0265 


G8 III-IV 


4819 


2.400 


-0.36 


2.718 


0.136 


0.091 


4 


6203 


HR 0296 


K0 III-IV 


4387 


1.750 


-0.48 


2.951 


0.139 


-0.181 


4 


6695 


HR 0328 


A3 V 


8390 


4.300 




1.531 


0.090 


0.071 


6 


7010 




K0 IV 








3.558 


0.283 


tt 


2 
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Table 1. Continued. 



HD 


Name 


Type" 


T c g b 


log/ 


[Fe/Hf 


D4000 


AD4ooo c 


Res. 


Runs 6 


7927 


HR 


0382 


F0 la 


7300 


0.400 




1.523 


0.090 


0.053 


6 


8424 


HR 


0398 


AO Vnn 


8455 


4.100 




1.400 


0.090 


-0.054 


6 


10307 


HR 


0483 


G2 V 


5898 


4.380 4 


0.04 4 


2.132 


0.086 


0.351 


5 


10380 


HR 


0489 


K3 III 


4133 


1.200 


-0.11 


4.323 


0.095 


0.561 


3 


10476 


HR 


0493 


Kl V 


5065 


4.500 1 


-0.20 5 


2.387 


0.137 


-0.053 


4 


10700 


HR 


0509 


G8 V 


5200 


4.560 


-0.37 


1.892 


0.086 


-0.276 


5 


10780 


HR 


0511 


K0 V 


5398 


4.600 1 


0.36 1 


2.319 


0.136 


-0.109 


4 


11004 






F7 V 


4842 






2.068 


0.088 


ft 
1 1 


2 6 


12929 


HR 


0617 


K2 III 


4480 


2.100 


-0.22 


3.106 


0.098 


-0.046 


6 


13043 






G2 V 


5730 


4.010 


0.07 


1.842 


0.086 


-0.078 


5 


13783 






G8 V 


5488 


4.500 


-0.55 


2.272 


0.136 


0.463 


4 


13974 


HR 


0660 


GO V 


5694 


4.600 


-0.33 s 


1.630 


0.091 


-0.140 


6 


14802 


HR 


0695 


Gl V 


5630 


3.920 


-0.17 


2.376 


0.094 


0.492 


6 


17378 


HR 


0825 


A5 la 


8530 


1.350 




1.480 


0.090 


0.237 


6 


17491 


HR 


0832 


M4 III 


3560 


0.600 




4.066 


0.105 


0.337 


6 


17709 


HR 


0843 


K5 III 


3905 


1.100 


-0.18 


4.033 


0.036 


-0.107 


6(7) 8 


18191 


HR 


0867 


M6 III 


3250 


0.300 




3.299 


0.099 


0.123 


6 


19373 


HR 


0937 


GO V 


5984 


4.000 2 


0.02 2 


1.872 


0.133 


0.154 


4 


19445 






G5 VI 


5830 


4.230 


-2.07 


1.400 


0.074 


0.123 


4 6 


19476 


HR 


0941 


K0 III 


4965 


2.550 


0.22 


2.826 


0.137 


-0.004 


4 


20041 


HR 


0964 


AO la 


9560 


2.130 




1.369 


0.090 


0.237 


6 


20630 


HR 


0996 


G5 V 


5528 


4.410 1 


0.10 1 


1.882 


0.085 


-0.232 


5 


20893 


HR 


1015 


K3 III 


4355 


1.950 


0.00 


3.128 


0.050 


-0.352 


1 4 6(2) 


22484 


HR 


1101 


F8 V 


6294 


4.010 1 


-0.04 4 


1.696 


0.062 


0.158 


5,6 


22879 






F9 V 


5780 


4.270 


-0.85 


2.180 


0.087 


0.629 


5 


23249 


HR 


1136 


K0+IV 


4820 


3.860 


0.02 


2.962 


0.097 


0.049 


6 


23439A 






Kl V 


5061 


4.500 1 


-1.02 1 


2.285 


0.135 


0.252 


4 


23439B 






K2 V 


4771 


4.500 1 


-1.05 1 


2.696 


0.137 


-0.047 


4 


25329 






Kl Vsb 


4862 


4.770 


— 1.69 4 


2.106 


0.087 


-0.138 


5 


26462 


HR 


1292 


F4 V 


6814 


4.120 


0.13 


1.463 


0.084 


0.015 


5 


26690 


HR 


1309 


F3 V 


6925 8 


4.210 s 


0.01 8 


1.521 


0.062 


0.067 


5,6 


26965 


HR 


1325 


Kl V 


5050 


4.520 


-0.03 


1.916 


0.204 


-0.647 


13 


27371 


HR 


1346 


K0 III 


4970 


2.600 


0.13 


2.614 


0.078 


-0.150 


4,6 


27697 


HR 


1373 


K0 III 


4983 


2.650 


0.13 


2.689 


0.078 


-0.060 


4 6 


28307 


HR 


1411 


K0 III 


4997 


2.690 


0.13 


2.473 


0.136 


-0.260 


4 


30455 






G2 V 


5671 


4.450 


-0.18 1 


1.783 


0.085 


-0.067 


5 


30649 






Gl V-VI 


5640 


4.320 


-0.37 6 


1.579 


0.085 


-0.208 


5 


30652 


HR 


1543 


F6 V 


6340 


4.290 


0.14 


1.583 


0.075 


0.037 


4 6 


34334 


HR 


1726 


K3 III 


4234 


2.150 


-0.16 


3.175 


0.140 


-0.395 


4 


34411 


HR 


1729 


GO V 


5848 


4.160 4 


0.10 4 


1.832 


0.034 


-0.008 


5 6(6) 


35620 


HR 


1805 


K3 III 


4156 


1.200 


0.21 


3.873 


0.143 


-0.013 


4 


36003 






K5 V 


4465 


4.610 


0.09 12 


3.179 


0.070 


0.063 


6(2) 


36162 


HR 


1832 


A3 Vn 


8260 


4.280 




1.427 


0.064 


-0.045 


6(2) 


37160 


HR 


1907 


OS TTTn 


4751 


2.900 


—0.55 


2.470 


0.095 


—0.123 


6 


38393 


HR 


1983 


F6 V 


6220 


4.170 


-0.07 


2.081 


0.092 


0.513 


6 


38751 


HR 


2002 


G8 III 


4751 


2.200 


0.03 


2.490 


0.060 


-0.468 


6,11(2) 


39587 


HR 


2047 


GO V 


5814 


4.400 1 


0.08 


1.671 


0.091 


-0.187 


6 


39970 


HR 


2074 


AO la 


9400 15 


1.430 




1.112 


0.059 


-0.034 


6,11(2) 


41117 


HR 


2135 


B2 lave 


17463 


2.700 




1.002 


0.055 


ttt 


6(2),11 


41636 


HR 


2153 


G9 III 


4690 


2.500 


-0.16 


2.777 


0.096 


-0.130 


6 


43318 


HR 


2233 


F6 V 


6069 


3.780 


-0.30 


1.756 


0.061 


0.165 


6,11(2) 


47205 


HR 


2429 


Kl IV 


4776 


3.350 


0.05 


3.495 


0.099 


0.522 


6 


47914 


HR 


2459 


K5 III 


3999 


1.400 


0.05 


4.097 


0.066 


0.067 


6,11(2) 


48433 


HR 


2478 


Kl III 


4465 


1.500 


-0.08 


3.134 


0.098 


-0.132 


6 
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Table 1. Continued. 



HD 


Name 


Type" 


T c g b 


logs 


[Fe/H] 6 


D4000 


AD4ooo c 


Res. 


Runs e 


48682 


HR 2483 


GO V 


5829 


4.080 


0.15 1 


1.717 


0.059 


-0.158 


6 11(2) 


49161 


HR 2503 


K4 III 


4276 


1.850 


0.23 


3.535 


0.100 


-0.206 


6 


49293 


HR 2506 


K0 III 


4635 


1.850 


-0.07 


2.990 


0.097 


-0.048 


6 


50778 


HR 2574 


K4 III 


4019 


1.350 


—0.13 


4.413 


0.107 


0.466 


6 


51440 


HR 2600 


K2 III 


4381 


2.200 


-0.35 


3.227 


0.061 


0.007 


6 11(2) 13 


54719 


HR 2697 


K2 III 


4381 


1.350 


0.24 


3.465 


0.099 


-0.142 


6 


54810 


HR 2701 


K0 III 


4697 


2.350 


-0.30 


2.650 


0.095 


-0.157 


6 


55575 


HR 2721 


GO V 


5640 


4.280 6 


-0.30 6 


1.628 


0.059 


-0.190 


6 11(21 


58207 


HR 2821 


K0 III 


4776 


2.350 


-0.10 


2.767 


0.096 


-0.073 


6 


58972 


HR 2854 


K3 III 


4023 


2.000 


—0.30 


3.534 


0.101 


—0.353 


6 


60522 


HR 2905 


M0 III-HIb 


3902 


1.200 


0.12 


3.996 


0.065 


-0.170 


6 11(21 13 


61935 


HR 2970 


K0 III 


4776 


2.200 


-0.02 


2.917 


0.096 


0.024 


6 


62721 


HR 3003 


K5 III 


3973 


1.450 


—0.01 


3.948 


0.075 


—0.104 


6 11(21 


64606 




G8 V 


5120 


4.110 


-0.92 


1.974 


0.060 


-0.017 


6 11(21 


66141 


HR 3145 


K2 III 


4264 


1.500 


-0.23 


3.435 


0.100 


-0.046 


6 


69267 


HR 3249 


K4 III 


4032 


1.600 


-0.36 


4.086 


0.104 


0.236 


6 


69897 


HR 3262 


F6 V 


6156 


4.150 


-0.41 4 


1.499 


0.090 


-0.037 


6 


70272 


HR 3275 


K5 III 


3931 


1.050 


-0.11 


4.116 


0.105 


0.015 


6 


72184 


HR 3360 


K2 III 


4669 


2.650 


0.34 


3.197 


0.092 


-0.097 


2 6 


72324 


HR 3369 


G9 III 


4864 


1.600 


0.17 


2.780 


0.053 


-0.137 


6 11(31 


73471 


HR 3418 


K2 III 


4480 


1.700 


0.09 


3.473 


0.062 


0.110 


6,11(21 


73593 


HR 3422 


G8 IV 


4717 


2.250 


-0.12 


2.695 


0.049 


-0.204 


6 11(41 


73665 


HR 3427 


K0 III 


4965 


2.350 


0.16 


2.762 


0.039 


-0.028 


1 8(21 11(21 


73710 


HR 3428 


K0 III 


4864 


2.100 


0.24 


2.853 


0.045 


-0.113 


6,8,11(31 


74377 




K3 V 


4913 


4.630 


-0.07 12 


2.757 


0.062 


0.022 


6,8 


74442 


HR 3461 


K0 III 


4651 


2.300 


—0.04 


2.997 


0.097 


—0.041 


6 


75732 


HR 3522 


G8 V 


5139 


3.700 


0.13 1 


2.608 


0.078 


0.051 


6 11(21 


79211 


BD+ 53 1321 


M0 V 


3769 


4.710 


-0.40 14 


3.036 


0.097 


0.223 


6 


82328 


HR 3775 


F6 IV 


6314 


4.010 


-0.10 4 


1.503 


0.059 


-0.022 


6 11(21 


83618 


HR 3845 


K3 III 


4240 


1.150 


-0.05 


3.813 


0.063 


0.190 


6,11(21 


84937 




F5 VI 


6230 


4.610 


-2.10 


1.265 


0.049 


-0.048 


6(21 11(21 


85503 


HR 3905 


K2 III 


4480 


2.150 


0.46 


3.428 


0.062 


-0.213 


6,11(2) 


86728 


HR 3951 


G4 V 


5675 3 


4.280 3 


0.11 3 


1.977 


0.092 


-0.008 


6 


88230 




K7 V 


3917 


4.680 


0.28 1 


2.961 


0.097 


0.038 


6 


88284 


HR 3994 


K0 III 


4971 


2.700 


0.22 


3.142 


0.061 


0.319 


6 11(21 


89010 


HR 4030 


G2 IV 


5670 


3.920 


-0.03 5 


1.849 


0.092 


-0.071 


6 


89449 


HR 4054 


F6 IV 


6334 


4.060 


-0.02 


1.482 


0.090 


-0.050 


6 


90508 


HR 4098 


Gl V 


5717 


4.370 


-0.23 1 


1.706 


0.060 


-0.091 


6 11(21 


94705 


HR 4267 


M5.5 III 


3335 


0.200 




3.323 


0.079 


-0.015 


6,11 


95272 


HR 4287 


K0 III 


4635 


2.150 


-0.06 


3.636 


0.100 


0.591 


6 


97907 


HR 4365 


K3 III 


4351 


1.600 7 


— 0.10 7 


2.736 


0.095 


—0.686 


6 


98230/1 HR 4374/5 


GO V 


5670 


4.410 


-0.21 1 


1.630 


0.091 


-0.207 


6 


101501 


HR 4496 


G8 V 


5360 


4.600 5 


-0.16 s 


1.999 


0.045 


-0.121 


3,6 


102328 


HR 4521 


K3 III 


4457 


2.350 


0.41 


3.669 


0.068 


0.038 


6,11(2) 


102870 


HR 4540 


F8 V 


6030 


4.050 


0.18 


1.602 


0.091 


-0.132 


6 


103095 


HR 4550 


G8 Vp 


5069 


4.600 4 


-1.32 4 


1.777 


0.043 


-0.122 


3,6 


106516 


HR 4657 


F5 V 


6140 


4.240 


-0.82 


1.492 


0.090 


0.021 


6 


108177 




F5 VI 


6050 


5.040 


-1.73 


1.342 


0.059 


0.008 


6,8 


109995 




AO V 


8300 


3.500 


-1.99 


1.380 


0.063 


-0.089 


6(2) 


110897 


HR 4845 


GO V 


5898 


4.290 6 


-0.44 6 


1.486 


0.031 


-0.141 


3(2) 


111721 




G6 V 


5103 


2.870 


-1.10 


2.010 


0.092 


0.001 


6 


113226 


HR 4932 


G9 II-III 


4994 


2.100 


0.08 


2.588 


0.051 


-0.117 


6(2), 11(3) 


113285 




M8 III 


2924 16 


0.000 




2.553 


0.061 


0.084 


6,8 


114710 


HR 4983 


GO V 


5900 


4.340 


0.15 


1.657 


0.037 


-0.163 


3,6(2) 



10 



J. Gorgas et al.: Empirical calibration of the A4000 A break 



Table 1. Continued. 



HD 


Name 


Type 11 


rri b 
J_ c ff 


logo 6 


[Fc/Hl 6 

1 1 J 


J — ^±UIJIJ 


Z^Dzinnn c 


Res. d 


Runs e 


114762 




F9 V 


5814 


4.260 


-0.75 4 


1.379 


0.042 


-0.190 


3 


114946 


HR 4995 


G6 V 


4930 


3.610 


-0.44 


2.768 


0.096 


0.239 


6 


114961 




M7 III 


3014 


0.000 


-0.81 


2.633 


0.095 


-0.046 


6 


115043 




Gl V 


6016 


4.510 


-0.08 1 


1.705 


0.046 


0.034 


3 


115617 


HR 5019 


G6 V 


5480 


4.520 


0.22 


1.976 


0.050 


-0.262 


3 


117176 


HR 5072 


G5 V 


5370 


3.800 


-0.21 


1.990 


0.033 


-0.091 


3,6(2),8(2) 


120136 


HR 5185 


F7 V 


6431 


4.270 


0.18 1 


1.497 


0.043 


-0.027 


3 


120452 


HR 5196 


K0.5 III-IIIb 


4760 


2.600 


0.03 


3.126 


0.098 


0.179 


6 


121146 


HR 5227 


K2 IV 


4428 


3.000 


-0.12 


3.638 


0.101 


0.345 


6 


121370 


HR 5235 


GO IV 


6036 


3.720 


0.28 6 


1.689 


0.091 


-0.069 


6 


122563 


HR 5270 


GO VI 


4626 


1.500 


-2.60 


1.632 


0.091 


-0.226 


6 


124897 


HR 5340 


K2 IIIp 


4276 


2.100 


-0.42 


3.166 


0.099 


-0.188 


6 


125454 


HR 5366 


G9 III 


4775 


2.500 


-0.10 


2.780 


0.042 


-0.061 


6,8(3) 


125560 


HR 5370 


K3 III 


4381 


1.650 


0.31 


3.454 


0.099 


-0.203 


6 


126327 




M7.5 III 


2943 11 


0.000 


-0.58 


2.402 


0.048 


-0.112 


6,8(2) 


129312 


HR 5480 


G8 III 


4834 


1.550 


0.00 


2.532 


0.094 


-0.304 


6 


131977 


HR 5568 


K4 V 


4493 


4.650 5 


0.01 5 


3.015 


0.067 


-0.103 


3 


132142 




Kl V 


5057 


4.500 5 


-0.55 


2.386 


0.049 


0.126 


3,11(2) 


132345 


HR 5582 


K3 III-IVp 


4374 


1.600 


0.42 


3.928 


0.102 


0.182 


6 


134083 


HR 5634 


F5 V 


6577 


4.380 


0.05 5 


1.438 


0.038 


-0.039 


3,6 


134439 




K0 V 


4990 


4.510 


-1.40 


2.291 


0.094 


0.282 


6 


135722 


HR 5681 


G8 III 


4834 


2.450 


-0.44 


2.553 


0.078 


-0.010 


4,6 


136028 


HR 5690 


K5 III 


4040 


1.900 


0.14 


4.301 


0.106 


0.298 


6 


136202 


HR 5694 


F8 III-IV 


6030 


3.890 


-0.07 


1.679 


0.091 


0.014 


6 


137391 


HR 5733 


F0 V 


7190 


4.140 


0.28 


1.450 


0.090 


-0.018 


6 


137471 


HR 5739 


Ml III 


3810 


1.100 


0.07 14 


3.708 


0.102 


-0.402 


6 


137759 


HR 5744 


K2 III 


4504 


2.050 


0.12 


3.714 


0.101 


0.363 


6 


139669 


HR 5826 


K5 III 


3915 


1.000 


-0.10 


4.504 


0.107 


0.375 


6 


140283 




F3 VI 


5650 


4.240 


-2.45 


1.202 


0.040 


-0.025 


3 


140573 


HR 5854 


K2 III 


4528 


2.100 


0.25 


3.404 


0.081 


-0.010 


4,6 


141144 


BD+01 3131 


K0 III 








3.006 


0.097 


tt 


6 


141680 


HR 5888 


G8 III 


4707 


1.950 


-0.30 


2.807 


0.091 


0.012 


2,6 


142091 


HR 5901 


K0 III-IV 


4751 


2.850 


0.02 


2.875 


0.079 


-0.077 


4,6 


142373 


HR 5914 


F9 V 


5718 


4.140 4 


-0.39 4 


1.478 


0.043 


-0.255 


3 


142860 


HR 5933 


F6 V 


6280 


4.270 


-0.25 


1.479 


0.090 


-0.039 


6 


142980 


HR 5940 


Kl IV 


4512 


2.550 


0.10 


3.625 


0.139 


0.300 


4 


143107 


HR 5947 


K3 III 


4318 


1.700 


-0.29 


3.450 


0.081 


0.093 


4,6 


143761 


HR 5968 


G2 V 


5835 


4.080 4 


-0.19 4 


1.662 


0.041 


-0.078 


3,6 


144872 




K3 V 


4739 


4.650 


-0.25 10 


2.880 


0.069 


-0.112 


3 


145148 


HR 6014 


K0 IV 


4849 


3.450 


0.13 


2.968 


0.065 


0.053 


4,6,11 


145328 


HR 6018 


K0 III 


4687 


2.500 


0.01 


2.790 


0.056 


-0.236 


4,6,11(2) 


145675 




K0 V 


5353 


4.490 1 


0.27 1 


2.777 


0.063 


0.365 


3 


147379A 


Gl 617 A 


M0 V 


3720 


5.000 2 


-1.50 2 


2.730 


0.059 


-0.040 


8(2) 


147379B 


Gl 617 B 


M3 V 


3247 


4.840 




2.023 


0.137 


-0.194 


8 


1 A *7RT7 


T 1 Tl £? 1 AO 

xlK 6103 


TVA TTT 

KU 111 


a q ao 
4893 


Z.35U 


a m 
U.Ol 


Z.t)20 


a i on 
0.139 


A 1 A O 

— U.148 


A 

4 


148513 


HR 6136 


K4 IIIp 


4046 


1.000 


0.20 


4.474 


0.144 


0.457 


4 


148783 


HR 6146 


M6 III 


3250 


0.200 


-0.06 


3.304 


0.053 


0.128 


6(2), 11(2) 


148816 




F9 V 


5815 


4.000 


-0.72 


1.534 


0.090 


-0.043 


6 


149161 


HR 6159 


K4 III 


3940 


1.450 


-0.13 


3.835 


0.147 


-0.248 


4 


149661 


HR 6171 


K0 V 


5133 


4.560 1 


0.01 1 


2.705 


0.061 


0.222 


3 


151203 


HR 6227 


M3 IHab 


3640 


0.700 


-0.10 14 


3.852 


0.103 


-0.005 


6 


152792 




GO V 


5647 


4.120 


-0.38 1 


1.728 


0.047 


-0.051 


3 


153210 


HR 6299 


K2 III 


4536 


2.250 


0.09 


3.123 


0.080 


-0.162 


4,6 


157089 




F9 V 


5749 


4.380 


-0.57 4 


1.557 


0.045 


-0.095 


3 
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Table 1. Continued. 



HD 


Name 


Type a 


T c ff b 




[Fe/H] 6 


D4000 


AD4ooo c 


Res. 


Runs e 


157214 


HR 6458 


Gl V 


5686 


4.240 6 


— 0.41 6 


1.655 


0.041 


-0.088 


3,6 


160693 




GO V 


5733 


4.220 


-0.67 5 


1.632 


0.085 


0.007 


5 


161797 


HR 6623 


G5 IV 


5390 


3.830 


0.23 


2.289 


0.093 


-0.053 


6 


161817 




A2 VI 


7450 


2.930 


-1.71 


1.377 


0.090 


—0.103 


6 


163506 


HR 6685 


F2 Ibe 


6400 


1.200 


-0.41 


1.425 


0.064 


-0.164 


6(2) 


164259 


HR 6710 


F3 V 


6750 


4.100 


-0.12 1 


1.428 


0.061 


-0.017 


5,6 


165195 




K3p 


4475 


1.200 


-2.70 


2.122 


0.083 


0.115 


2 5 


165760 


HR 6770 


G8 III-IV 


4957 


2.000 


-0.05 


2.689 


0.030 


0.026 


4 5(2) 6(7) 7 12 


165908 


HR 6775 


F7 V 


5883 


4.230 


-0.50 6 


1.481 


0.090 


-0.134 


6 


166620 


HR 6806 


K2 V 


4941 


4.650 


-0.25 1 


2.560 


0.065 


—0.041 


5,6 


167006 


HR 6815 


M3 III 


3640 


0.700 


-0.17 14 


4.157 


0.106 


0.300 


6 


167042 


HR 6817 


Kl III 


4949 


3.050 


-0.06 


2.804 


0.106 


0.136 


2(3),4 


168720 


HR 6868 


Ml III 


3810 


1.100 


o.oo 14 


4.086 


0.104 


—0.024 


6 


168775 


HR 6872 


K2 III 


4520 


1.700 


0.07 


3.094 


0.142 


-0.199 


4 


170153 


HR 6927 


F7 V 


5968 


4.360 


-0.33 1 


1.521 


0.062 


-0.106 


5,6 


172958 


HR 7030 


B8 V 


11300 


3.750 




1.162 


0.059 


-0.091 


6 11(2) 


175638 


HR 7141 


A5 V 


8150 


3.900 




1.561 


0.090 


0.077 


6 


175743 


HR 7148 


Kl III 


4635 


2.450 


-0.09 


3.304 


0.093 


0.279 


2,6 


175751 


HR 7149 


K2 III 


4697 


2.200 


0.03 


3.035 


0.142 


0.008 


4 


175865 


HR 7157 


M5 III 


3759 11 


0.500 


0.14 14 


3.461 


0.100 


-0.575 


6 


176301 


HR 7171 


B7 III-IV 


13100 


3.500 




1.129 


0.090 


-0.041 


6 


176411 


HR 7176 


K2 III 


4687 


2.100 


0.17 


2.880 


0.069 


-0.261 


4 11(2) 


180928 


HR 7317 


K4 III 


3969 


1.300 


-0.35 


4.641 


0.103 


0.655 


2 6 


181984 


HR 7352 


K3 III 


4443 


2.400 


0.39 


4.053 


0.098 


0.419 


5 


182572 


HR 7373 


G8 IV H<51 


5555 


3.380 


0.38 5 


2.571 


0.095 


0.297 


6 


184406 


HR 7429 


K3 III 


4428 


2.450 


0.22 


3.095 


0.142 


—0.433 


4 


185144 


HR 7462 


K0 V 


5133 


4.500 5 


-0.24 5 


2.468 


0.088 


0.142 


5 


185859 


HR 7482 


B0.5 Iae 


21860 


2.800 




1.264 


0.059 


ttt 
1 1 1 


6 11(2) 


186408 


HR 7503 


G2 V 


5831 


4.060 


0.12 5 


1.761 


0.085 


-0.100 


5 


186427 


HR 7504 


G5 V 


5639 


4.060 


0.07 5 


2.161 


0.061 


0.165 


5(2) 


187691 


HR 7560 


F8 V 


5968 


3.990 


0.10 6 


1.546 


0.084 


-0.207 


5 


187923 


HR 7569 


GO V 


5665 


4.210 


0.06 1 


1.998 


0.086 


0.029 


5 


188056 


HR 7576 


K3 III 


4355 


2.150 


0.42 


3.705 


0.063 


-0.065 


3,6 


188512 


HR 7602 


G8 IV 


5062 


3.700 


-0.35 


2.516 


0.095 


0.157 


6 


188727 


HR 7609 


G5 lb var 


5685 


1.670 


0.00 


1.826 


0.060 


-0.224 


6 11(2) 


190360 


HR 7670 


G6 IV + M6 V 


5248 


3.890 


0.26 


2.330 


0.094 


-0.194 


6 


190406 


HR 7672 


Gl V 


5989 s 


4.350 s 


-0.07 s 


1.849 


0.086 


0.160 


5 


190603 


HR 7678 


B1.5 Iae 


18800 


2.410 




1.495 


0.090 


ttt 


6 


195593 


HR 7847 


F5 lab 


6600 


1.950 


0.12 


2.169 


0.093 


0.648 


6 


195633 




GO V 


5830 


3.780 


-1.12 


1.842 


0.085 


0.376 


5 


197076 


HR 7914 


G5 V 


5760 


4.240 


-0.08 


1.944 


0.086 


0.113 


5 


198478 


HR 7977 


B3 Iae 


16450 9 


2.100 9 


-0.23 9 


1.123 


0.089 


ttt 

1 1 1 


6 


199478 


HR 8020 


B8 Iae 


11200 


1.900 




1.186 


0.089 


ttt 
1 1 1 


6 


1 99580 




KO TTT-TV 


5039 


3.500 


—0.10 


2.778 


0.042 


0.243 


1 9(9) 3 6 


201626 




G9p 


4941 7 


2.000 7 


-1.50 7 


2.218 


0.055 


0.249 


3 


201891 




F8 V-VI 


5890 


4.400 


-1.23 


1.658 


0.085 


0.228 


5 


203344 


HR 8165 


K0 III-IV 


4669 


2.600 


-0.09 


2.931 


0.042 


-0.049 


1(2),2,3 


205153 




GO IV 


5910 


3.700 


-0.06 


1.784 


0.078 


0.045 


8 


205650 




F6 V 


5810 


4.450 


-1.24 


1.505 


0.078 


0.063 


8 


207260 


HR 8334 


A2 la 


9450 


2.090 




1.582 


0.090 


0.440 


6 


208906 




F8 V-VI 


5930 


4.280 


-0.82 


1.463 


0.078 


-0.052 


8 


210027 


HR 8430 


F5 V 


6413 


4.200 


o.oo 1 


1.602 


0.132 


0.089 


4 


213470 




A3 la 


8700 15 


1.380 




1.345 


0.090 


0.124 


6 


215648 


HR 8665 


F7 V 


6010 


3.830 


-0.40 


1.623 


0.132 


0.032 


4 
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Table 1. Continued. 



HD 


Name 


Type" 


rj-i b 
± c ^j- 


logo 6 


[Fe/Hf 






Res. d 


Runs 6 


217476 


HR 8752 


G4 la 


5100 


1.500 7 


0.07 1 


1.635 


0.091 


-0.949 


6 


219734 


HR 8860 


M2 III 


3730 


0.900 


0.27 14 


3.990 


0.104 


-0.003 


6 


222368 


HR 8969 


F7 V 


5959 


3.940 


-0.25 6 


1.644 


0.133 


-0.009 


4 


232078 




K4-5 III-II 


4000 


0.400 


-1.60 


3.989 


0.100 


0.384 


6,11(2) 




\ji III 


1X0 V 


3 /D» 


4. (00 


— 0.33 


2.959 


0.091 


0.14b 


r 





BD+ 17 4708 


F8 VI 


5900 


3.930 


-1.90 


1.368 


0.078 


0.064 


8 




BD+ 19 5116A 


M4 V 


3200 


4.910 




2.310 


0.070 


0.158 


6(2) 




BD+ 19 5116B 


M6 V 


2950 


5.060 




2.238 


0.098 


0.438 


6 




BD+ 43 0044B 


M6 V 


3107 


5.080 




2.013 


0.101 


-0.007 


6 




BD+ 44 2051A 


M2 V 


3544 


4.850 




2.360 


0.049 


-0.232 


6,8,11(2) 




BD+ 56 1458 


K7 V 


4070 


4.700 


-0.18 13 


3.098 


0.063 


0.089 


6,8 




Gl 699 


M5 V 


3153 


5.000 




1.992 


0.066 


-0.093 


6(2) 




Luyton 789-6 


M7e 


2747 


5.090 




1.360 


0.078 


-0.182 


8 




Ross 248 


M6e 


2799 


5.120 




1.832 


0.146 


0.228 


6 



a Sources for spectral types are, in order of preference: W94, G93, the Bright Star Catalog (Hoffleit 1982) and 
the Simbad database. For the cluster stars (except for Coma and Hyades) we list positions in the HR diagram 
(SGB: subgiant branch; GB: giant branch; HB: horizontal branch, red clump stars; AGB: asymptotic giant 
branch) . 

b Atmospheric parameters from the electronic table of W94, with the following additional sources and modifi- 
cations: (1) Cayrel de Strobel et al. (1997). (2) Alonso et al. (1996). (3) Gratton et al. (1996). (4) Mean from 
sources 1, 2 and 3. (5) Mean from sources 1 and 2. (6) Mean from sources 1 and 3. (7) Thevenin (1998). (8) 
Marsakov and Shevelev (1995). (9) Gies and Lambert (1992). (10) Zakhozhaj and Shaparenko (1996). (11) 
Dyck et al. (1996). (12) Computed from the indices Fe5270 and Fe5335 using the fitting functions of G93 (in 
terms of (V — K)). (13) Computed from the indices Fe5270 and Fc5335 using the fitting functions of W94 (in 
terms of T e ff)- (14) Computed from the index Fe5270 using the fitting functions of W94 (in terms of T e g). (15) 
From spectral type using table III of Johnson (1966). (16) Extrapolation in the relation T a g versus spectral 
type of Ridgway et al (1980) — these temperatures differ from those given by W94 — . 
c Total random errors in the A4000 A break measurements. 

d Residuals from the derived fitting functions (observed — predicted). Stars not used to derive the empirical 
functions are marked with | (non cluster members, see G93), ft (lacking required atmospheric parameters) 
or f|| (emission-line supergiants) . 

e Run number(s) in which the star was observed. Number of repeated observations within each run is giving 
in parenthesis. 



sources of the stellar parameters and their associated er- 
rors are given in the original papers G93 and W94. In this 
work we assume that these errors are random and, thus, 
their effect in the fitting procedure is minimized through 
the use of a library containing a large number of stars. The 
other type of errors are those associated to D4000 measure- 
ments, which are the subject of this section. Undoubtly, 
an accurate knowledge of the errors is essential to guaran- 
tee the validity of the final product of this work, i.e., the 
fitting functions of the break. 



Since, apart from the cluster members, most of the 
stars of the Lick/IDS library are bright, and considering 
the low signal-to-noise ratio required to measure the D4000 
with acceptable accuracy, systematic errors are the main 
source of uncertainty. 



5.1. Random errors 



(i) Photon statistics and read-out noise. With the aim 
of tracing the propagation of photon statistics and read- 
out noise, we followed a parallel reduction of data and 
error frames. For a detailed description on the estima- 
tion of random errors in the measurement of line-strength 
indices we refer the interested reader to Cardiel et al. 
(1998). Starting with the analysis of the photon statis- 
tics and read-out noise, the reduction package RePuiE is 
able to generate error frames from the beginning of the 
reduction procedure, and properly propagates the errors 
throughout the reduction process. In this way, important 
reduction steps such as flatfielding, geometrical distortion 
corrections, wavelength calibration and sky subtraction, 
are taken into account. At the end of the reduction pro- 
cess, each data spectrum S(Xi) has its associated error 
spectrum a(Xi), which can be employed to derive accu- 
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Run 
IXU.I1 


Dates 


r TYil Anrid" 


spectrograph 


Detector 


uisp. 


AA c 


r vv nivi 


1 


1-3 Nov 1991 


WHT 4.2m 


ISIS blue arm 


EEV#6 


2.70 


3500- 


-6850 


16.0 


2 


8-10 Jun 1993 


CAHA 3.5m 


CTS 


TEK#12 


3.46 


o/OU- 




7.2 


3 


2-6 Jul 1994 


INT 2.5m 


IDS 235mm 


TEK#3 


3.30 


3700- 


-7100 


5.2 


4 


9-10 Aug 1994 


CAHA 3.5m 


CTS 


TEK#12 


3.46 


3700- 


-7240 


8.1 


5 


10-14 Sep 1994 


INT 2.5m 


IDS 235mm 


TEK#3 


3.30 


3700- 


-7100 


5.8 


6 


13-19 Feb 1995 


JKT 1.5m 


RBS 


TEK#4 


0.91 


3600- 


-4530 


2.8 


7 


15 Feb 1995 


INT 2.5m 


IDS 235mm 


TEK#3 


0.85 


3634- 


-4502 


2.1 


8 


17-19 Dec 1995 


WHT 4.2m 


ISIS blue arm 


TEK#1 


2.90 


3740- 


-6700 


12.3 


9 


15-16 Jan 1996 


CAHA 3.5m 


CTS 


TEK#11 


0.88 


3654- 


-4554 


2.5 


10 


4 Aug 1996 


WHT 4.2m 


ISIS blue arm 


TEK#2 


1.53 


3654- 


-5219 


4.9 


11 


28 Mar-1 Apr 1995 


SPM 2.12m 


B&Ch 


TEK 


4.00 


3500- 


-7592 


8.3 


12 


3 Mar 1995 


INT 2.5m 


IDS 500mm 


TEK#3 


0.75 


3670- 


-4440 


4.2 


13 


17-18 Nov 1996 


CAHA 3.5m 


CTS 


SITe#6a 


1.10 


3569- 


-5765 


3.4 



Q JKT, INT and WHT (La Palma, Spain), CAHA (Calar Alto, Spain), SPM (San Pedro Martir, Mexico) 

6 dispersion (A/pixel) 

c wavelength range (A) 

d measured spectral resolution (A) 



rate index errors. The errors in the break are computed 
by (Cardiel et al. 1998) 



A 2 [D 



4000jphoton 



T r G 2 Tb + TbOjFr 



with 



and 



N p 



[^(Ai)], 



(4) 



(5) 



(6) 



where the subscripts b and r correspond, respectively, to 
the blue and red bandpasses of the break (p refers indis- 
tinctly to b or r), S(Xi) and cr(A;) are the signal and the 
error in the pixel with central wavelength Ai, O is the dis- 
persion (in A/pixel) assuming a linear wavelength scale, 
and N p is the number of pixels covered by the p band (in 
general, fractions of pixels must be considered at the bor- 
ders of the bandpasses). We have checked that the above 
analytical formulae exhibit an excellent agreement with 
numerical simulations. For the whole sample, the error of 
a typical observation introduced by these sources of noise 
is (A_[D 4 ooo] phot?n ) = 0.038. _ 

(ii) Flux calibration. During each run we observed a 
number (typically around 5) of different spectrophotomet- 
ric standard stars (from Massey et al. 1988 and Oke 1990). 
The break was measured using the average flux calibration 
curve, and we estimated the random error in flux calibra- 
tion as the r.m.s. scatter among the different D4000 values 



obtained with each standard. The typical error introduced 
by this uncertainty is (A [D 40 oo]fi ux ) = 0.034. 

(iii) Wavelength calibration and radial velocity correc- 
tion. These two reduction steps are potential sources of 
random errors in the wavelength scale of the reduced spec- 
tra. Radial velocities for field stars were obtained from the 
Hipparcos Input Catalogue (Turon et al. 1992), which in 
the worst cases are given with mean probable errors of ~ 
5 km s" 1 (~ 0.07 A at A4000 A). For the cluster stars, we 
used either published radial velocities for individual stars, 
if available, or averaged cluster radial velocities (Hesser 
et al. 1986: M3, M5, M10, M13, M71, M92, NGC6171; 
Friel 1989: NGC188; Friel & Janes 1993: M67, NGC7789; 
Turon et al. 1992: Coma, Hyades). Typical radial velocity 
errors for the cluster stars are ^ 15 km s" 1 (~ 0.2 A at 
A4000 A) . To have an estimate of the random error intro- 
duced by the combined effect of wavelength calibration 
and radial velocity, we cross-correlated fully calibrated 
spectra corresponding to stars of similar spectral types. 
The resulting typical error is 20 km s _1 , being always be- 
low 75 km s _1 . This translates into a negligible error of 
(A [D 4 ooo] wavolcngth > = 0-003. However, it may be useful to 
estimate the importance of this effect when measuring the 
break in galaxies with large radial velocity uncertainties. 
As a reference, using the 18 spectra displayed in Fig. ||, a 
velocity shift of ~ 100 km s _1 translates into relative D4000 
errors always below 1%. Furthermore, for K0 III stars we 
obtain A [D4000] wavelength — 1-56 x 10~ 4 Aw, where Av is 
the velocity error in km s" 1 (this relation only holds for 
Aw < 150 km s _1 ; for Aw in the range from 150-1000 
km s _1 the error increases slower, and remains below 0.1). 

(iv) Additional sources of random errors. Expected 
random errors for each star can be computed by adding 
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quadratically the random errors derived from the three 
sources previously discussed, i.e., 



A [D4000] expected — 



A [D4000] photon + A [D 40 oo]flux + A 2 [D 40 oo] wavelength ■ (7) 



However, additional (and unknown) sources of ran- 
dom error may still be present in the data. Following 
the method described in Gonzalez (1993), we compared, 
within each run, the standard deviation of the D4000 mea- 
surements of stars with multiple observations with the ex- 
pected error A[D4ooo]cxpcctcd- For those runs in which the 
standard deviation was significantly larger than the ex- 
pected error (using the -F-test of variances with a signifi- 
cance level a = 0.3), a residual random error A res iduai was 
derived and added to all the individual stellar random 
errors: 



A 2 [D4ooo]random — 

A [D4000] expected + A [D4000] residual- (8) 

It is worth noting that this additional error was only 
needed for some runs. In the particular case of run 6, with 
a large number (54) of stars with multiple observations, 
the agreement between expected and measured error was 
perfect. 

5.2. Systematic effects 

The main sources of systematic effects in the measure- 
ment of spectral indices in stars are spectral resolution, 
sky subtraction and flux calibration. 

(i) Spectral resolution. We have examined the effect of 
instrumental broadening in the break by convolving the 
18 spectra of Fig. || with a broadening function of variable 
width. The result of this study indicates that, as expected, 
the break is quite insensitive to spectral resolution. As a 
reference, for a spectral resolution of 30A (FWHM) the 
effect in the break is below 1%. Therefore, given the res- 
olutions used in this work (last column in Table |^) no 
corrections are needed in any case. 

(ii) Sky subtraction. Since the field giant and dwarf 
stars of the library are bright, the exposure times were 
short enough to neglect the effect of an anomalous sub- 
traction of the sky level. However, most of the cluster 
stars are not bright, being necessary exposures times of 
up to 1800 seconds for the faintest objects. In addition, 
the observation of these stars, specially those in globular 
clusters, were performed with the unavoidable presence of 
several stars inside the spectrograph slit, which compli- 
cated the determination of the sky regions. In Cardiel et 
al. (1995) we already studied the systematic variations on 
the D4000 measured in the outer parts of a galaxy (where 
light levels are only a few per cent of the sky signal) due 
to the over- or under-estimation of the sky level. We refer 
the interested reader to that paper for details. Although 
there is not a simple recipe to detect this type of system- 
atic effect, unexpectedly high D 4 ooo values in faint cluster 
stars can arise from an anomalous sky subtraction. 

(Hi) Flux calibration. Due to the large number of runs 
needed to complete the whole library, important system- 
atic errors can arise due to possible differences among 
the spectrophotometric system of each run. In order to 
guarantee that the whole dataset is in the same system, 
we compared the measurements of the stars in common 
among different runs. Since run 6 was the observing run 
with the largest number of stars (including numerous mul- 
tiple observations) and with reliable random errors (see 
above), we selected it as our spectrophotometric reference 
system. Therefore, for each run we computed a mean offset 
with run 6, which was introduced when it was significantly 
different from (using a t test). It is important to high- 
light that differences between a true spectrophotometric 
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system and that adopted in this work may still be present. 
Therefore, we encourage the readers interested in the pre- 
dictions of the present fitting functions, to include in their 
observations a number of template stars from the library 
to ensure a proper correction of the data. 

5.3. Final errors 

The comparison of measurements of the same stars in dif- 
ferent runs also provides a powerful method to refine the 
random errors derived in Eq. ^| We followed an iterative 
method which consistently provided the relative offsets 
and a set of extra residual errors to account for the ob- 
served scatter among runs (see Cardiel 1999 for details). 

As mentioned before, the data sample was enlarged 
by including 43 stellar spectra from the Pickles' (1998) 
library. The random errors in the D4000 indices measured 
in this subsample were estimated from the residual vari- 
ance of a least-square fit to a straight line using all the 
stars (except supergiants) with T c g > 8400 K (they follow 
a tight linear relation in the D4ooo-0 plane). The typical 
error in Pickles' spectra was found to be 0.036. 

6. D4000 measurements 

Table |l| (columns 7 and 8) lists the final D4000 measure- 
ments and the associated random errors. Some sample 
spectra, exhibiting a diversity in spectral types, metal- 
licities and gravities, are displayed in Fig. ||, whereas in 
Fig. |I] we show the break behaviour with effective tem- 
perature for the whole sample. It is clear from these plots 
that temperature is the main parameter governing D4ooo- 
The effect of metallicity is clearly noticeable in panel ||(b) 
and by the position of globular cluster stars in Fig. |J. Also, 
some gravity dependences are also observed, especially be- 
tween hot dwarf-giants and supergiants, and between cold 
dwarfs and giants. 

The definition of the D4000 given in Eq. |l| resembles 
that of a color. However, the peculiar combination of v 
and A translates into the introduction of a wavelength 
weighting of the flux (Eq. ||). In order to facilitate the 
computation of the break, we have studied the effect of a 
redefinition of the index, namely 



B 



4250 (1+z) 



fx dA 



4000 



_ ./ 4050 (1+z) 



3950 (1+z) 



(9) 



fx dA 



3750 (1 + z) 



Comparing this and the previous definition of D4000 
(Eq. ||) for all the stars in run 6, we find that no single star 
deviates more than a 1% from the theoretical predicted 
ratio 



D 



4000 



B 



4000 



= 1.1619, 



(10) 



obtained for a constant f\. Therefore the above ratio can 
be safely used to convert between both break definitions. 

7. The fitting functions 

The main aim of this work is to derive empirical fitting 
functions for the D4000 in terms of the stellar atmospheric 
parameters: effective temperature, metallicity and surface 
gravity. After some experimentation, we decided to use 
9 = 5040/T e ff as the temperature indicator, being [Fe/H] 
and log g the parameters for the metallicity and gravity. 
Following the previous works of G93 and W94, the fitting 
functions are expressed as polynomials in the atmospheric 
parameters, using two different functional forms: 



D4000 (9, [Fe/H] , log 9) =p(0, [Fe/H] , log g) 



and 



D 



4000 \ 



?, [Fe/H], log g) = const. 



p(e,[Fc/H]4ogg) 



(11) 



(12) 



where p is a polynomial with terms up to the third order, 
including all possible cross-terms among the parameters: 
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p(9, [Fe/H], log g) =]Tc fc 0* [Fe/H] j (logs)' 



(13) 



fc=0 



with < i + j + I < 3. 

The polynomial coefficients were determined from a 
least squares fit where all the stars were weighted accord- 
ing to the D4000 observational errors listed in Table 1. Note 
that this is an improvement over the procedure employed 
by G93 and W94 for the Lick indices. 

Obviously, not all possible terms are necessary. The 
strategy followed to determine the final fitting function is 
the successive inclusion of terms, starting with the lower 
powers. At each step, the term which yielded a lower new 
residual variance was tested. The significance of this new 
term, as well as those of all the previously included coeffi- 
cients, was computed using a t-test (i.e. from the error in 
the coefficient, we tested whether it was significantly dif- 
ferent from zero). Note that this is equivalent to perform- 
ing a F-test to check whether the unbiased residual vari- 
ance is significantly reduced with the inclusion of the addi- 
tional term. Following this procedure, and using typically 
a significant level of a = 0.10, only statistically signifi- 
cant terms were retained. The problem is well constrained 
and, usually, after the inclusion of a few terms, a final 
residual variance is asymptotically reached and the higher 
order terms are not statistically significant. Throughout 
this fitting procedure we also kept an eye on the residuals 
to assure that no systematic behaviour for any group of 
stars (specially stars from any given cluster or metallicity 
range) was apparent. 

After a set of trial fits, it was clear that tempera- 
ture is the main parameter governing the break. Unfor- 
tunately, the behaviour of the D4000 could not be repro- 
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Fig. 3. Sample spectra of stars observed in run 6. Effective temperatures are given in parenthesis. Panel (a) is a 
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duced by a unique polynomial function in the whole tem- 
perature range spanned by the library, forcing us to di- 
vide the temperature interval into several regimes. The 
derived composite fitting function is shown in Fig. ^. In 
Table || we list the corresponding coefficients and errors, 
together with the typical error of the N stars used in each 
interval (of yp = N/ Y^i=i a i 2 )i the unbiased residual vari- 
ance around the fit (c 2 td ) and the determination coefficient 
(r 2 ). 

In the high temperature regime (9 < 0.75, T e g > 6700 
K) a dichotomic behaviour for dwarfs and giants on one 
side, and supergiants on the other, is clearly apparent. 
Therefore we derived different fitting functions for each 
gravity range. For the first group the amplitude of the 
break is quite constant and only the linear term in 9 is 
statistically significant (note that we subdivide this range 
in two intervals to achieve a better fit). The independence 
on mctallicity is naturally expected (see section 2) but 
note that an important fraction of the stars in this range 
either lack of a [Fe/H] estimation or are restricted to the 
solar value. 

The behaviour of the cool stars (0.75 < 9 < 1.3, 
3900K < T cff < 6700 K) is more complex and [Fe/H] 



terms are clearly needed. On the other hand, no grav- 
ity term is significant. However, whilst for the giant stars 
D4000 increases with 9 all the way up to 9 ks 1.3, for higher 
gravities it reaches a maximum at 9 w 1.1 and then lev- 
els off. Furthermore, separate fits for dwarfs and giants in 
this T e ff range (with a gravity cutoff around 3 — 3.5) yield 
residual variances that are significantly smaller than the 
variance from a single fit. Hence, we have derived differ- 
ent fitting functions for dwarfs and giants. This dichotomic 
behaviour of the break is not surprising since its strength 
is quite dependent on the depth of the CN bands (Fig. 1) 
which also shows a similar behaviour (G93) due to the on- 
set of the dredge-up processes at the bottom of the giant 
branch. In Fig. || we show in detail the fitting functions 
derived for each gravity group in this temperature range. 

Concerning the cold stars (9 > 1.3, T e3 < 3900 K), 
the difference between giants and dwarfs is quite evident 
and two fitting functions have been derived (see also sec- 
tion 2). Again, the metallicity terms are not significant, 
although this may be, at least in part, due to the paucity 
of input metallicities in this range. It must be noted that 
the different fitting functions have been constructed with 
the constrain of allowing for a smooth transition in the 
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predicted D4000 indices among the different T e g and grav- 
ity ranges. 

In Fig. fj] we plot the residuals from the fits as a func- 
tion of effective temperature, metallicity and gravity. Note 
that no trends are apparent with any of these parameters. 
We have also checked for systematic residuals within any 
of the star clusters. Except for an unexplained negative 
offset for the Coma stars (AD4000 = 0.09, not due to an er- 
ror in the adopted metallicity) , no systematic offsets have 
been found. For the 420 stars used in the fit, we derive an 
unbiased residual standard deviation er st d = 0.160. This 
must be compared with the typical error in the D40007 
Ctyp = 0.064. Therefore, the residuals are, in the mean, 
a 2.5 factor larger than what should be expected solely 
from measurement errors. Since we are quite confident 
that this latter errors are realistic (see Section |5|), and 
although some scatter may arise from the fact than the 
fitting functions are not able to reproduce completely the 
complex behaviour of the D4000, most of the extra scatter 



must arise form uncertainties in the input atmospheric pa- 
rameters. For example, the residual D4000 scatter of 0.248 
for the cool giants (at 9 = 1.0 and [Fe/H] = 0.0) can be 
fully explained by the combined effect of a 166 K uncer- 
tainty in T e g and a 0.29 dex error in [Fe/H], both con- 
sistent with the typical errors found by Soubiran, Katz 
& Cayrel (1998) when comparing atmospheric parameters 
from the literature. Another quantitative measurement of 
the quality of the present fitting functions is the deter- 
mination coefficient for the whole sample r 2 = 0.96. This 
indicates that a 96% of the original variation of the break 
in the sample is explained by the derived fitting functions. 

Since the goal of this work is to predict reliable D4000 
indices for any given combination of input atmospheric pa- 
rameters, we have investigated, using the covariance ma- 
trices of the fits, the random errors in such predictions. 
These errors are given in Table |] for some representative 
sets of input parameters. Note that, as it should be ex- 
pected, the uncertainties are smaller for near-solar metal- 
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Table 3. Parameters of the empirical fitting functions in 
each temperature and gravity range 



Hot stars 


0.12 < 


9 < 0.63 


3 < logg < 4.5 


polynomial fit 






N = 39 


c : 


0.6548 


± 0.0222 


(Jtyp = 0.041 


9 : 


1.340 


± 0.049 


a std = 0.046 








r 2 = 0.96 



Warm stars 0.63 < 9 < 0.75 3 < logg < 4.5 



polynomial fit 

c : 1.823 ± 0.140 
9 : -0.5068 ± 0.2030 


N = 23 

CT t yp = 0.051 

<r std = 0.034 
r 2 = 0.35 


Hot supergiants 0.19 < 9 < 0.69 < logg < 3 


polynomial fit 

c : 0.8613 ± 0.0328 
9 3 : 1.849 ± 0.211 


N = 18 

(Jtyp = 0.041 

<T std = 0.079 
r 2 = 0.87 


Cool dwarfs 0.75 < 9 < 1.08 3 < logg < 5.1 


exponential fit 

CO 

9 

[Fe/H] 
9 [Fe/H] 
9 2 

9 2 [Fe/H] 


(const. = 0.9) 

-8.154 ± 2.007 

13.45 ± 4.28 
-17.06 ± 7.28 

38.12 ± 15.80 
-4.769 ± 2.269 
-20.69 ± 8.52 


N = 161 

(Jtyp = 0.068 
cr st d = 0.178 
r 2 = 0.88 


Cool giants 0.75 < 9 < 1.30 < logg < 3.5 


exponential fit 

CO 

9 

[Fe/H] 
9 [Fe/H] 
9 2 

9 2 [Fe/H] 


(const. = 0.9) 

-5.665 ± 1.563 

9.279 ± 2.766 
-3.273 ± 2.274 

7.322 ± 4.266 
-3.080 ± 1.218 
-3.694 ± 1.986 


N = 176 

(Jtyp = 0.081 

o-std = 0.261 
r 2 = 0.85 


Cold dwarfs 1.08 < 9 < 1.83 4.5 < logg < 5.2 


exponential fit 

CO 

9 
9 2 


(const. = 1.15) 

-2.908 ± 2.949 
6.535 ± 4.243 
-2.976 ± 1.506 


N = 15 

(Jtyp = 0.078 

cr std = 0.194 
r 2 = 0.83 


Cold giants 1.30 < 9 < 1.72 < logg < 1.2 


polynomial fit 

c : 9.525 ± 0.742 
9 : -4.094 ± 0.486 


N = 15 

(Jtyp = 0.083 
cr std = 0.211 
r 2 = 0.92 



licities. Interestingly, although the library does not include 
a high number of 0-B stars, the predicted indices at the 
hot end of the star sample are rather reliable. 

8. Summary 

We have derived a set of empirical fitting functions de- 
scribing the behaviour of the break at A4000 A in terms 
of the atmospheric stellar parameters: effective tempera- 
ture, metallicity and surface gravity. This calibration can 




Fig. 6. Details of the fitting functions in the mid- 
temperature range for a) giant (log g < 3.5) and b) dwarf 
(logg > 3) stars. See caption to Fig |{| Error bars for the 
D4000 measurements are also shown 



be easily incorporated into stellar population models to 
provide accurate predictions of the D4000 for composite 
systems. In a forthcoming paper we will analyze the mea- 
surements of the break in old stellar populations at the 
light of the predictions of such models. Considering the 
volume covered by the employed stellar library (Lick/IDS 
+ Pickles' hot subsample) in the stellar parameter space, 
the derived fitting functions suit the requirements to pro- 
vide accurate D4000 predictions for populations with ages 
larger than about 0.1 Gyr, and — 1 < [Fe/H] < +0.5 dex 
(see also W94). Note however that, since the break could 
be contaminated by nebular emission in galaxies with on- 
going star formation, it is necessary to include this effect 
in the population modelling in order to predict reliable 
D4000 indices for such stellar populations. It should also 
be noted that the applicability of the derived fitting func- 
tions is safe as far as the abundance ratios of the library 
stars reflect those in the modeled stellar populations. 

In order to facilitate the usage of the present D4000 
fitting functions we have written a FORTRAN subroutine, 
available from: 

frittp : //www . ucm . es/inf o/Astrof /D4000/D4000 . html . 

This routine computes the value of the D4000 as a function 
of the input stellar parameters T e g, [Fe/H] and logy. The 
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Fig. 7. Residuals of the derived fitting functions (observed 
minus predicted) against the three input stellar parame- 
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of the error bar is twice the unbiased residual standard 
deviation. 



code performs smooth interpolations among temperature 
and gravity ranges, providing also an estimate of the error 
in the predicted index, 
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Table 4. Uncertainties in the predicted D 40 oo values for 
some indicative sets of stellar parameters. For the cool 
stars (T c ff < 6000K) the term giants refers also to super- 
giant stars. 





[Fe/H] 


dwarfs-giants 


supergiants 


30000 




0.015 


0.032 


1 9nnn 

1ZUUU 




U.UUo 


U.UZo 


7500 




0.010 


0.044 








cn 3 trf"t; 


6000 


0.5 


0.045 


0.156 


6000 


0.0 


0.024 


0.099 


6000 


-1.0 


0.031 


0.075 


6000 


-2.0 


0.043 


0.105 


5000 


0.5 


0.103 


0.079 


5000 


0.0 


0.044 


0.043 


5000 


-1.0 


0.082 


0.053 


5000 


-2.0 


0.119 


0.073 


4000 


0.5 




0.167 


4000 


0.0 


0.125 


0.088 


4000 


-1.0 




0.173 


4000 


-2.0 




0.322 


3200 




0.094 


0.082 
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